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Abstract 

For a bilinear form obtained by adding a Dirac mass to a positive definite moment 
functional in several variables, explicit formulas of orthogonal polynomials are de- 
rived from the orthogonal polynomials associated with the moment functional. Ex- 
plicit formula for the reproducing kernel is also derived and used to establish certain 
■ inequalities for classical orthogonal polynomials. 
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1 Introduction 



Let U d denote the space of polynomials in ^-variables, and let u be a moment 
functional, denoted by (u,p) for p e U d , for which orthogonal polynomials 
exist. We define a new functional v by adding a Dirac mass to u, 

(v,p) = (u,p) + \p(c), ceR d , AgM, peU d , (1) 



1 Partially supported by Ministerio de Ciencia y Tecnologfa (MCYT) of Spain and 
by the European Regional Development Fund (ERDF) through the grant MTM 
2005-08648-C02-02, and Junta de Andaluci'a, Grupo de Investigacion FQM 0229. 
The work of the fourth author is supported in part by NSF Grant DMS-0604056. 



Preprint submitted to Elsevier 



2 February 2008 



and study orthogonal polynomials with respect to the functional v. 

In the case of one variable, this problem first arose from the work of A. M. 
Krall ([5]) when he studied the orthogonal polynomials that are eigenf unctions 
of a fourth order differential operator considered by H. L. Krall ( [BUT] ), and 
showed that the polynomials are orthogonal with respect to a measure that 
is obtained from a continuous measure on an interval by adding masses at 
the end points of the interval. In [3], Koornwinder studied the case that the 
measure is a Jacobi weight function together with additional mass points at 
1 and —1; he constructed explicit orthogonal polynomials and studied their 
properties. Uvarov ([ID]) considered the problem of orthogonal polynomials 
with respect to a measure obtained by adding a finite discrete part to another 
measure; his main result expresses the polynomials orthogonal with respect to 
the new measure in terms of the polynomials orthogonal with respect to the 
old one. More generally, one can consider perturbations of quasi definite linear 
functionals via the addition of Dirac delta functionals, orthogonal polynomials 
in such a general setting has been studied extensively in recent years (see, for 
instance [8], and the references therein). 



The purpose of the present paper is to study this problem in several variables. 
After a brief section on notations and preliminaries in the next section, we 
state and prove our main results in Sections 3. The result gives a necessary and 
sufficient condition for the existence of orthogonal polynomials with respect to 
the linear functional v defined by ([T]), and expresses orthogonal polynomials 
with respect to v in terms of the orthogonal polynomials with respect to u. 
Furthermore, we can also express the reproducing kernel of polynomials with 
respect to v in terms of the kernel with respect to u. Our formula on the 
reproducing kernel implies an inequality on the orthogonal polynomials with 
respect to u. Even in the case of one variable, it leads to new inequalities 
on classical orthogonal polynomials, which are stated in Section 4. Finally, in 
Section 5, we consider the example of orthogonal polynomials on the unit ball. 



2 Orthogonal polynomials in several variables 

In this section we recall necessary notations and definitions about orthogonal 
polynomials of several variables, following [T]. 

Throughout this paper, we will use the usual multi-index notation. For a = 
(cti, . . . , a d ) e Nq and x = (xi, . . . , x d ) G M d , we write x a = x" 1 ■ • • x a d d . 
The integer |a| = «i + • • • + ay is called the total degree of x a . The linear 
combinations of x a , \a\ = n, is a homogeneous polynomial of degree n. We 
denote by V% the space of homogeneous polynomials of degree n in d variables, 
by the space of polynomials of total degree not greater than n. It is well 
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known that 

dim IE? = ( "~ r "| a nd dim 7^ 





Let {/i a } Qe Nd be a multi-sequence of real numbers, and let u be a real valued 
functional defined on V by means of 

(u,x a ) = u{x a ) = fj, a , 

and extended by linearity. Then, u is called the moment functional determined 
by {/Ua,} QeN d. If (u,p 2 ) > 0, Wp G LT^ and p ^ 0, then the moment u is 
called positive definite and it induces an inner product accordingly by (p, q) : = 
(u,pq), Vp, g G n d . A typical example of a positive moment functional is 
an integral with respect to a positive measure dfi with all moments finite, 
(u,p) = J R dp(x)dfi. 

A polynomial P G LT^ is called an orthogonal polynomial with respect to 
u if (u,PQ) = for all Q G n^_ 1 . Let denote the space of orthogonal 
polynomials with respect to the u. We are interested in the case when u admits 
a basis of orthogonal polynomials; that is, dim V^f = for all n. This happens 
whenever u is positive definite. In general, we call a moment functional quasi 
definite if it admits a basis of orthogonal polynomials. 

For the study of orthogonal polynomials it is often convenient to adopt a 
vector notation ( [HfTTj ). Let {P2}\a\=n denote a basis of V*- Let the elements 
of {a G N d : \a\ = n} be ordered by a±, a 2 , ■ ■ ■ , a. r d according to a fixed 
monomial order, say the lexicographical order. We then write the basis of 
as a column vector 

— (P a )\a\=n = (Pan Pa 2 i ■ ■ ■ i Pa d ) ■ 

n 

Using this notation, the orthogonality of P™ can be expressed as 

(u P P T ) = H 5 

where H n is a matrix of size x r~|. That u is quasi definite is equivalent 
to that H n is invertible for all n > 0. If u is positive definite, then we can 
choose a basis so that H n is the identity matrix for all n; in other words, we 
can choose the basis to be orthonormal. 

Let u be a quasi definite moment functional and let {Pnjn^o be a sequence of 
orthogonal polynomials with respect to u. The reproducing kernel of V£ and 
LT^ are denoted by P n (-, •) and K n (-, •), respectively, and are given by 

n 

P k (x,y) = Fl(x)H k 1 ¥ k (y) and K n (x, y) = £ P k (x, y). 

k=0 
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These kernels satisfy the usual reproducing property; for example, 

p(x) = (u,K n (x,-)p(-)) = (u,p(-)K n (-,x)), p G U d , 

which shows, in particular, that these functions are independent of a particular 
choice of the bases. The kernel K. n (x,y) satisfies an analog of the Christoffel- 
Darboux formula, and plays an important role in the study of orthogonal 
Fourier expansions. For this and further properties of orthogonal polynomials 
of several variables, see [T]. 



3 Krall— type orthogonal polynomials in several variables 

Let u be a quasi definite moment functional defined on U d . We define a new 
moment functional v as the perturbation of u given by 

(v,p) = (u,p) + Ap(c), Vp{x) G n d , 

where A is a non zero real number and c G M. d is a given point. Our first result 
gives a necessary and sufficient condition for v to be quasi definite. 

Theorem 1 The moment functional v is quasi definite if and only if 

A n := 1 + AK n (c,c) ^ 0, n > 0. 

Furthermore, when v is quasi definite, a sequence of orthogonal polynomials, 
{Qn}n>o, with respect to v is given by 

Q„0) = ¥ n (x) - K n _!(c, x) P n (c) n > 0, (2) 

where {P n }n>o denote orthogonal polynomials with respect to u and K_!(-, •) : = 
0. 

PROOF. First we assume that v is quasi definite and {Q n } n >o is a sequence 
of orthogonal polynomials with respect to v. Since u is quasi definite, the 
leading coefficient of P n is an invertible matrix. Hence, by multiplying Q n by 
an invertible matrix, if necessary, we can assume that Q n — P n G V n -\ for 
n > 0. This shows, in particular, that Qo = Po- Furthermore, since {P n }n>o is 
a basis of U d , for each n > 1 we can express Q n in terms of P n . Thus, there 
exist constant matrices M" of size r d x rf such that 

n-1 

Q n (x)=F n (x) + Y,MJ l F, l (x), 

i=0 
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where, by the orthogonality of P n and the definition of v, 

Ml = (u,Q n Fj)H^ = [(v,Q n Fj) - A Q n (c) Pf (c)] Hr\ 
for < % < n — 1. Since (v, Q n Pf ) = for i < n — 1, we conclude then 

n-l 

Q n (x) = P n (x) - A Q n (c) X! p f(c) fff 1 ^0) = " A Q„(c) K n _!(c, x). 

i=0 

Evaluating the above expression at x — c, we obtain 

Q n (c)[l + AK n _ 1 (c,c)]=P n (c). (3) 

Recall that A& = 1 + AKjt(c,c). If A n _i = for some value of n > 1, then 
P n (c) = by © and, furthermore, A n = A n _i + AP^(c) i/- 1 P n (c) = 0. Thus, 
we conclude that P n+ i(c) = which, however, contradict to the fact that P n 
and P n +i cannot have a common zero (PQ, p. 113). This shows that A„ ^ 0, 
Vn > 0, and also that (j2J) holds whenever X n ^ 0. 

Conversely, if all A n are non-zero, then we can define the polynomials Q n by 
§Z§ and the above proof shows then that Q n is orthogonal with respect to v. 
Since Q n and P n has the same leading coefficient matrix, it is evident that 
{Qn}n>o contains a basis of U d . □ 

Remark 2 If u is a quasi definite moment functional, then v is quasi definite 
except when A belongs to an infinite discrete set of values. If u is positive 
definite and we choose A > 0, then it is easy to see that v is also positive 
definite (see (j3J below). 

Let Q n be as in the theorem, we define H k = (v, Qj[) and denote by 

n 

P k (x, y) := QT(x) H, 1 Q k (y) and K n (x, y) := J2 ?k(x, y), 

k=0 

the reproducing kernels associated with the linear functional v. We will derive 
an explicit formula of these kernels when u is positive definite. 

Whenever u is positive definite, we choose {P„}„>o in the Theorem [T] as a 
sequence of orthonormal basis with respect to u, so that H n = (u,F n Fl) = 7 Pb , 
the identity matrix of size r n . 

Proposition 3 Let u be a positive definite linear functional. Then, fork > 0, 
H k = I rk + -^F k (c)F k (c) T and H k l = I rjt - ^F k (c) P fc (c) T . (4) 

A k 



PROOF. Since we choose P„ so that H n = I rn , it follows from (El) that 
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H k = (v,Q k QZ) = (v,Q k Fl) = (u,Q k ¥l) + XQ k (c)F k (c) T 

= I rk + -^F k (c)F k (c) T . 

Assuming that the inverse is of the form H* 1 = I r . k - 6F k (c) P fc (c) T , and 
using the fact that ¥ k (c) T F k (c) = P n (c, c), a quick computation shows that 
HkH^ 1 = I Tk is equivalent to A — \ k ~i5 — X5P n (c, c) = 0. Using the fact that 
P n (c, c) = K n (c, c) — K n _x(c, c), it is easy to see that 5 = X/\ k . 

Theorem 4 Let u be a positive definite linear functional. Then, for k > 0, 

P k (x, y) = P k (x, y)--^- K k (x, c) K fe (c, y) + — — K fe _i(x, c) K fc _ x (c, y). (5) 
Furthermore, for n > 0, 

Kn(x, y) = Kn(x, y) - — K n (x, c) K n (c, y). (6) 

PROOF. Let 7 fc := A/A fe . By © and (0}, it follows readily that 

P*(z, 2/) = p fc(^ I/) - 7fe [Pfc(x, c)K fc _i(ac, c) + P fc (y, c)K fc _i(c, x) 

+P k (x, c)P k (y, c)] + 7fe7fc-iPfc(c, c)K fc _i(x, c)K fc _i(y, c) 
= P k (x, y) - 7fcK fc (x, c)K fc (j/, c) 

+ Tfet 1 + 7fc__iP fc (c, c)]Kjfc_i(a:, c)K fc _i(y, c), 

as the first square bracket is equal to Kfe(x, c)Kfc(y, c) — K^_i(x, c)K&_i(y, c). 
Since the definition of 7^ leads readily to 1 + 7fe_iPfe(c, c) = 7^-1, this proves 
(JSJ). Summing over dSJ) for = 0, 1, . . . , n proves (J6]). □ 

We note that, in the case of one variable, the formula ([6]) has appeared in [2], 
whereas the formula (jSJ) appears to be new even in one variable. 



4 An application of formula (J5J) 

If u is positive definite and A > 0, then v is also positive definite. As a result, 
H n , hence H~ x are positive definite matrices by (j3J) for n > 0. In particular, 
P n (x, x) is nonnegative. In fact, if d > 2 and the linear functional u is centrally 
symmetric (see [I] for definition), then P n (x,x) is strictly positive for all x, 
except for n odd and x = 0. As a consequence, we see that implies that 

P n (x, x) - A [K„(x, c)] 2 + [K n _x(x, c)] 2 > 

An A n -1 
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for all x G M. d and for all A > 0. Recall that X n = 1 + AK n (c, c). Taking the 
limit A — > oo, we obtain an inequality which we state as a proposition. 



If x — c, then the two sides of (jTJ) are equal. This inequality tends out to 
be non-trivial even in the case of one variable. Notice that in one variable, 
P n (x,x) = \p n (x)] 2 , where p n is the orthonormal polynomial. Let us specify 
the inequality (j7j) in the cases of classical orthogonal polynomials of Jacobi 
and Laguerre. We use the standard notation P^ x, ^\ a,j3 > —1, for Jacobi 
polynomials, which are orthogonal with respect to (1 —x) a (l + x)@ on [—1, 1], 
and a > —1, for the Laguerre polynomials, which are orthogonal with 
respect to x a e~ x on [0, oo). Using the formulas in [9], especially (4.3.3), (4.5.3) 
and (5.1.1), (5.1.13), the inequality (0) for the Jacobi polynomials with c = 
1, and for the Laguerre polynomials with c = 0, respectively, becomes the 
following: 

Proposition 6 For n > 1 and x G [—1, 1], 



[Pj a ' 0) (x)] 2 n + (3 [P^X 1,P \x)) 2 > n + a + P + 1 [P^\x)] 2 

Pt ,P \l) + 2n + a + f3 + l P^+ 1,/3) (l) ~ 2n + a + (3 + 1 pt +m (l) 



As far as we are aware, these inequalities are new. For example, in the case 
of Chebyshev polynomials or a = /? = - 1/2 in the Jacobi polynomials, the 
inequality becomes 

2 „ 1 /sin(n- ±)6\ 2 1 /sin(n + i)fl\ 2 

2 cos 2 nd + — i > — i ^- , < 9 < tt. 

2n - 1 \ sin | J ~ 2n + l\ sin f J ~ ~ 

5 An example: Krall— type orthogonal polynomials in the unit ball 

Let B d denote the unit ball of M. d . We consider the inner product 



Proposition 5 For n > 1, and x, c G M. d , 



P n (x,x) + 



[K^foc)] 2 > [K n (x,c)] 2 
K n _ x (c,c) ~ K n (c, c) 



(7) 



For n > 1 and x G [0, oo) 



[L^{x)f [Lt\ l \x)} 2 y [4" +1) (*)] 2 
Lt\0) 4- + i 1} (0) ~ 4 a+1) (0) 




7 



where fi > —1/2, and c M = + ^p)/(7r d / 2 r(yU + |)) is the normalization 
constant so that (1, 1) M = 1. As an example for our general results, we add 
the mass point at the origin, and consider the inner product 

(/,ff) = (/,^ + A/(0M0) ) A>0. (8) 

We now use (T3]) in Theorem [T] to find an orthogonal basis for (•,•). 

Let Ti^ denote the space of spherical harmonic polynomials of degree n in d 
variables. Let Y™, 1 < v < dim7Y^, be an orthonormal basis for Ti^ in the 
following. An orthonormal basis for (/, g)^ is given explicitly by ([H p. 39]) 



P U X ) 



h n 



h ' n 2j+ ^\2\\x\\ 2 -l)Yr 2j (x), 0<j<n/2, (9) 



where pj ' 2 denote the orthonormal Jacobi polynomials, and hj v is 

r 1 2 

the normalizing constant given by kj v = (d / 2) n _ 2 j / (fJ> + ^f-)n-2j, in which 
(a)^ := a(a + 1) . . . (a + k — 1) denotes the shifted factorial. Since FJ 1-2 - 7 is a 
homogeneous polynomial, Y™~ 2: >(0) = unless its degree is zero, that is, unless 
n = 2j. Consequently, -Pj^(O) = unless j = n/2 and n is even. Notice that 

when n = 2j, W = 1. Hence, it follows that PR u (0) = p„ 2 ' 2 (—1) if n is 
even and Pj 1 y (0) = in all other cases. Consequently, if we define polynomials 

Q].„ by 

n n ( \ - \ P lA x ) -PnK n _i(x,0), if n is even 
I j> W 5 otherwise 

where p n := Ap^'^^-l)/^ + AK n _x(0, 0)), then according to © in The- 
orem [H : 1 < v < dim^_ 2j -,0 < 2j < n} constitutes an orthogonal 
basis with respect to the inner product ([8]). 

To make the expression for transparent, we assume [i > and make use of the 
following explicit formula for the reproducing kernel K n (-, •) in [12J, 

K n (x,y) = A% £ P^'^- 1 \x -y + y/l- \\x\\*y/l - \\y\\H)(l - t^dt, 

where \i > (see [12] for the case /x = 0) and 

2T(n + i)r(/i + ^)r(n + 2/i + d) 



7rV2r( M )r(n + /x + f )r(2/i + d + 1) ' 



We set y = in this formula and follow through a sequence of manipulations 
of formulas. First we use P, (4.5.3)] to write p^ At+2 ' /i+ 2 1 ^ ag a gum Q f 
Gegenbauer polynomials, so that we can apply [TJ Theorem 1.5.6] to get ride of 
the integral (for the terms of odd degree Gegenbauer polynomials, the integrals 
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are automatically zero), the result is a sum of |_§ J terms of Jacobi polynomials 
upon using the first formula on [TJ, p. 27], which we can use [9l (4.5.3)] again 
to sum up. The final result is the following identity, 



K n (x,0) 



2 J p\ 2 >r" 2 



(l-2||x|| 2 ), 



where [x\ denotes the integer part of x. In particular, setting x = gives 



11^ 



K„,(0,0) 



(A*+|)LiJ 



2 J 2 >A* 2 



(1) 



(12) 



Substituting these formulas into (fTUj) gives a basis of explicit orthogonal poly- 
nomials with respect to the inner product in (jSJ). Furthermore, using (fTTj) 
and (1121) in the formula ([H]), we obtain a compact formula for the reproduc- 
ing kernel K n (-, •) associated with (•, •) in (|S]). We sum up these results as a 
proposition. 

Proposition 7 For the inner product (•, •) in i/ie polynomials Q n j V , 1 < 
< dim?^_ 2 - ; < j < n, in fflUT) /brm an orthogonal basis of degree n. 
Furthermore, the reproducing kernel ofU^ with respect to (■, ■) is given by 



K n (x, y) = K n (x, y) - d n P^ i 
where the constant d n is given by 

A 



LfJ 



1-2 1/ 



AK„(0,0) 



+ |)Lf J 
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